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Abstract. 


This paper presents a class of new exf)licit second order accurate finite 
difference schemes for the conputaiton of weak solutions of hyperbolic conservation 
laws. '<hese highly nonlinear schemes are obtained by applying a nonoscillatory 
first order accurate scheme to an appropriately modified flux function. The so 
derived second order accurate schemes achieve high resolution while preserving 
the robustness of the original nonoscillatory first order accurate scheme. 

Numerical experiments are presented to demonstrate the perfomance of these 


new schemes. 


- 1 - 


1. Introduction . 

In this paper we consider numerical approximtions to weak solutions of the 
initial value problem CIVP) for hyperbolic systems of conservation laws. 

uCx,0) = i^(x), -«* < X < ». Cl.l) 

Here u(x,t) is a cloumn vector of m unknowns, and f(u), the flux, is a 
vector valued function of m conponents. Cl.l) is called hyperbolic if all 
eigenvalues a^Cu),..., SjjjCu) of the Jacobian natrix A(u) 

ACu) = f^ a. 2a) 

are real and the set of right eignevectors F?-(u),..., Ff'^Cu) is complete. 

I’fe assume that the eigenvalues {a^ (u)> are arranged in a nondecreasing order 

cp-(u) ^ a^Cu) S ... i a"*Cu). (1.2b) 


Vfe consider systems of conservation laws (1.1) that possess an entropy 
functicHi U(u), defined as follows: 

(i) U is a convex function of u, i.e., > 0, 

(ii) U satisfies 


Vu 



U.3a] 


where F is some other function called entropy flux. Admissible weak solutions 
of (1.1) satisfy, in the weak sense, the following inequality: 
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U(u)^ + FCu)^ ^ 0 a. 3b) 

(see [11]). The inequality (1.3b) is called an entropy condition. 

In the following we shall discuss numerical approximations to weak solutions 
of Q.l) v^ch are obtained by (2k+ l)-point explicit schemes in conservation 
farm 



where 


^4 




n 

'j+k‘ 


jVV^u) 


(1.4b) 


Here = vCjAx, nAt), an ? is a numerical flux function. We require the 
numerical flux function to be consistent with the flux f(u) in the following 
sense: 

?(u,...,u) = f(u). (1.4c) 


Vfe say that the difference scheme (1.4) is consistent with the entropy condition 
(1.3b) if an. inequality of the following kind is satisfied: 


U? - X 



a.5a) 


vhere Uj = U(Vj), ■ ^^'^j-k+1** * * ’'^'j+k^ ’ here ? is a numerical entropy 

flux, consistent with the entropy flux F(u), i.e. 
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f(u,...,u) = FCu). CL- 5b) 

Vfe trim now to discuss the question of convergence of the finite difference 
solution of Q.4) to weak solutions of the conservation laws (1.1). Since the 
finite-difference scheme is nonlinear and the computed solutions are certainly 
not smooth, therefore L 2 - stability of a consistent finite-dlference scheme 
does not iJiply convergence. One can establish convergence of finite - difference 
solutions of Q.4) to weak solutions of (1.1) when the following conditions 
are satisfi^; 

(i) The total variation with respect to x of the finite-difference solutions 
is uniformly bounded with respect to t, At and Ax. 

(ii) The finite-difference scheme (1.4) is consistent with the entropy ocndition 
(1.3b) for all entropy functions of Q.l). 

-Ciii) The entropy condition (1.3b) Implies uniqueness of the solution to the IVP 

( 1 . 1 ). 


Using compactness arguments one can deduce from condition (i) the existence 
of convergent subsequences. The conservation form (1.4) and condition (ii) 
imply that each limit solution is a weak solution which satisfies the entropy 
condition Q.3b). When the entropy condition implies uniqueness of the IVP 
(conditon (iii)) then all subsequences have the same limit solution, and con- 
sequently the finite-difference scheme is convergent, (see [2] , [9] , [lOJ ) . 


It seems possible to satisfy conditions (i) and (ii) by adding a hefty 
amount of artificial viscosity to the finite difference scheme (1.4). The addi' 
tional viscosity terms damp possible oscillations in the computed solution, and 
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HBke the convergence process siinuilate the zero-dlssipatlon-'lijnit which is used 
to select the unique physicedly relevemt weak solution. Unfortunately, viscosity 
represents an irretrievable loss of infomation and therefore the addition of 
artificial viscosity brings about some deterioration in resolution. 

In this paper we describe a new method to design finite-difference schemes 
that satisfy conditions Ci) and (ii), but scxi second-order accurate and have 
high resolution. 


2 
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. yonotonicit7/ in the scalar case. 

In this section we consider the IVP for a scalar conservation law. 

+ fCu)^ = + a(u)u^ = 0 , a(u) = ^ C2.1a) 

u(x,0) = ♦Cx), — < X < «», C2.1b) 

where (x) is assumed to be of bounded total variation. Every 
weak solution of the scalar IVP (2.1) which satisfies the entropy 
condition has the following mono tonicity propert y as a function of t 

i) No new local extrema in x nay be created. 

ii) The value of a local minirraim is nondecreasing, the value of a local mximum 
is nonincreasing. 

It follows from this monotonicity property that the total variation in 
X, T/(u(t)), of u(x,t), is nonincreasing in t, i.e., 


Ti/(u(t2)) TVCu(t^)), for all t2 ^ t^. 


( 2 . 2 ) 


We consider now eocplicit (2k + 1) point finite difference schemes in con- 
serv’ation form (1.4) apprcximting (2.1) 


n+1 n n n ^ n n ..n . 

Vj = H( '^j-k+1’ * • • ’''j+k^ " j ^^''j-k+1’ • * * 


(2.3a) 


15 #1 \ T 


and denote (2.3a) in an operator form as 
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= L-v^ 


C2.3b) 


We say that the finite difference scheme (2.3) is total variation nonincreas- 
ing CIV!rn if for al_l V of bounded total variation 


TV(L-v) ^ TVCv) 
v/here 

TVCU> = j ; 

•j :_oo J X/ <- 

here, and throughout this paper, we use the standard notation 


(2.4a) 


(2.4b) 




“ ^ 3 ’ 


(2.5) 


V7e say that the fiiiite-difference scheme (2.3) is monotonicity preserving 
if ‘thie finite difference operator L is monotonicity preserving, i.e., if v 
is a monotone mesh f miction so is L*v. 


We say that the finite difference scheme (2.3) is a monotone scheme if 
H in (2.3a) is a monotone nor decreasing function of each of its 2k + 1 
arguments. 

The following theorem states the hierarchy of these properties. 

T heoi~em 2.1. (i) A monotone scheme is TVNI. (ii) A TVNI scheme is monotonicity 

preser\'ing. 
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Proof : (i) It was proven by B. Keyfitz in [8] that monotone schemes form an 
f, ^-contractive semigroup, i.e. 


j|L*v - L-z!| $ i|v - z|l 

for all i, -sirniTBble v and z; Here | |uj | 

1 '1 

imrrediately from applying (2.0a) to v and z 


(2.6a) 


00 

= ^ lu.j. (2.4) follows 

• _ j 

3=-« 

= T*v, (i.e. z- = V..T for all i). 
’ J D+1 


(ii) Let (2.3) be a TVIIE scheme and let v be a monotone mesh function of 
bounded total variation, and denote w = L*v. Since L la.s a finite support 
of 2k + 1 points it is sufficient to prove that w is lacfiotone for all v of 
the form 


constant = 

monotone 

constant 


j 

i j ^ 5 • 

j >. J+ 


TV(v) = |vp^ ■" \^1 


C2.6b) 


Vfe prove (ii) by negation. Suppose w is not monotone, then it has at 
least one local miniinum and one local, maximum. Denote by and the values 
of the first two successive local extrema, then 

TV(w) 5 |Vj^ - -v^l > TV(v), 

which contradicts the assurrption that the scheme is TVNI. This conpletes the 
proof of Tneorem 2.1 
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Monotone scheires approximte solutions of the viscous modified equation 

'^y.^x * ^ ^ (2.7a) 

p(u,X) = I i^H.(u,u,...,u) - X^a^Cu)] (2.7b) 

2X^U=-k J 

p(u,X) ^ 0, eCu,X) ^ 0 (2.7c) 


to second or'der accuracy; Since g(u X) $ 0, 
first order accurate; in C2.7c) denotes 


monotone schemes ar>e necessarily 

(see [6 


]) 


Since monotone scliemes are T\/NI, there exist convergent subsequences for all 
initial data of bounded total variation. Each limit is a weak solution of C2.1) 
that satisfies Oleinik’s entropy cond5.tion (see C 8 ]). Since Oleinik’s entropy 
condition uiplies uniqueness of the IVP (2.1), we conclude that all subsequences 
converge to the same limit, and therefore the scheme is convergent (see [2 3). 

Let us consider new the scalar constant coefficient case a(u) = constant 
in (2.1). A linear finite-difference approximtion 


.n+1 


k 

= l 


t=-k 




Cj^ = const. 


C2.8a) 


is monotonicity preserving if and only if 


> 


c, i 0 


-k $ 1 ^ k. 
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(see [4 3), Hence any l inear monotonicity preserving scheme, and therefore any 
TVI«I linear scheme, is a monotone scheme, and consequently first order accurate. 

We renark that the previous statement does not exclude the possibility of 
having nonlinear monotonicity preserving and TVNI schemes that are second order 
accurate (and consequently are not monotone schemes). In fact the schemes presented 
in £7 3 and [6 3 are monotonicity preserving (at least in the constant coefficient 
case) and second order accurate. 


It is the purpose of this paper to present new high-resolution second order 
accurate TVNI scheires. These new schemes are generated by converting known 
3-point first order accurate TVI'C schemes into new 5-point second order accurate 
TVNI schemes. Both the 3-point schemes and the new 5-point schemes can be re- 
written in the form 


v^‘^^ = L-v^ 


C2.9a) 


(L-v)j = Vj + “ ^-,j-l/2^j-l/2'' 


(2.9b) 


where A. .1 is defined in (2.5) and 
^T 


^+,j+l/2 ■ ‘^-,1-1/2 ' (2.9o) 


The following Lenira states conditions on thie coefficients (2.9c) which are 


sufficient to ensoure that the scheme (2.9 ) is TVNI. 
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coefficients C In (2.9c) satisfy the Inequalities 
®-,j+l/2 - ° • °+,J+l/2 - ® (2.10a) 

®-,j+l/2 * ^+,1+1/2 - ^ • (2.10b) 

then the scheme (2.9) Is TVNI. 

Proof,: Denote w • L.v and subtract (2.9b) at j - 1 from (2.9b) at j « 1+1 

to obtain 

^1+1/2'' “ ®-,i-l/2^i-l/2^ ^^"^-.1+1/2"^+. 1+1/2^ ^1+1/2^ ^+,1+3/2 ^i+3/2''* 

( 2 . 11 ) 


Ey (2.10) all the coefficients in C2.ll) are non negative, therefore 


^^i>l/2'^* ^ (1 - ”^+^+l/2^'^i+l/2'^ * ^-,i-l/2^^i-l/2^* 


^+,i+3/2'^i+3/2^l 


Sunniing C2.12) for -» < i < « we obtain 


( 2 . 12 ) 


TVCw) i J ^ J " ^-,i+l/2 ” ^+,i+l/2^^^i+l/2''^ 


i=“ 


1=-' 


C. i_i/2lVi/2vl + .1 C+,i+3/2lW;''l = J JW''! 

-t-i •- J./ iz«eo ' X-“** 


wbi< 2 h shows that (2.4) is satisfied. The equality 5s obtained by changing the 
suitmtion index in the last two sums in the RHS of t:i inequality. 



- M - 


In the next section we shall use Leinrra 2.2 to design second order accurate 
T/NI schenres. We reirark ttet any 3 point finite difference scheme in conservation 
form with a differentiable numerical flux can be rewritten as (2.9), in the follow- 
ing way: It follows from the mean value theorem that there exist and C_ 

such that 


xl?(vj.v.,i) - rcv.,vj)] = -C/v.,Vj,3^) 

- fCVjjVj)] = - 


(2.13a) 

(2.13b) 


Expressing the numerical flux values in C2.3a) with k = 1 by (2.13) results 
in the form C2.9). 


3. Second order accurate TVllI schemes. 


Let us consider a general 3-point finite difference scheme in conservation 
form (1.4) with a numerical flux f of the form 


r(Vj,Vj+i> = A.J v] 


where 


•’ 2 


[fCv.+l) - fCvj)]/Aj+l V 


aCvj) 


when A. .1 V 4 0 

3+7 


when A. .1 V = 0 


(3. la) 


C3.Ib) 


Here Q(x) is some function, which is often referred to as the coefficient of 
numerical viscosity. 


Lenma 3.1 Let Q(x) in (3.1a) satisfy the inequalities 


|x| ^ Q(x) $ i 


for 0 ^ lx| $ y s 1, 


(3.2) 


then the finite-difference scheme (1.4) with (3.1) is TVNI under the CFL-like 
restriction. 


Xmax-|a^‘ ll $ y 
3 3+j 


(3.3) 


Proof: Using the notation 
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where a.J. is (3. lb), v;e rewrite (3.1a) as 
3+7 


xr.^l = xf(v.,v.^3^) = XfCv^) - ^ 


(3.4b) 


and similarly 


XF. 1 = Xf(v. - ,v.) = XfCv.) - yfv.JL + Q(^iJi)l • 

3“2 3 3 *^L 3-7 3 2 J 3 2 


C3.4c) 


Substituting (3.4) for numerical flux values in (1.4) we get the fom (2.9) 


= V? - «r.^ - ?..i) = vV t ^ 3,1 


n 


V 


(3.5a) 


where 


|[Q(v..l) . V..1J ...1 


n - n , 

V = V. + 


^4v"-C_ 1. IV. 


,n 


3 -+'3-2 3-7 -'3-2" 3^2 


Since 


C..54 - C.,.4 = QC -v.,| 


(3.5b) 


(3.5c) 


it follows from C3.2) and (3.5) that conditions (2.10) of Lemma 2.2 are satisfied 
under the CFL restriction C3.3) and therefore the finite-difference scheme (3.1) 
is TVNI. 


The second order accurate Lax— V7endro.ff scheme has the numerical flux (3.1) 
2 

wT-th Q(x) = X , i.e. , 
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XW _ 1 

^j+1/2 ^ 2 


[fCVj) + - 


1 

X 


^''j+l/2^^^j+l/2^^* 


C3.6ay 


Clearly a numerical flux of a second order accurate scheme 
to satisfy 


^j+1/2 


^+ 1/2 


= O(A^), 


C3.6b) 


for all smooth solutions of (1.1); here A is the discretization paraineter. 
When Q(x) is constrained by C3.2), then the 3-point scheme (3.1) is only 
first order accurate, for 


^j+1/2 


- ?LW 
■ j+1/2 




C- (v. 


j+1/2 




and therefore 


(3.7) 




j+1/2 " ^j+1/2 




2\ 


j+1/2' 


" ^'^j+1/2^ ^‘l^j+i/2 


v^‘l = 0(A). 


V7e describe new how to convert a 3-point first cir'.er accurate TVNI scheme to 

a 5-point second order accurate TTNI scheme. Consider the application of a 

<-M 

3-point first order accurate TVIiE scheme C3.1) to modiried mesh values 
of the original flux fCu): Set 


- f (v^) + 7 . gj ” g<''i-i.'i’ 'i+i) <3. 8a) 

where g is a function which will be specified below. Tie modified numerical 
flux f J+1/2 * ^^^j-l»^j»^j+l»^j+2^ obtained by rerlacing f(v^), f(v^^^) 
in (3.1a), (3.1b), by the modified values: 
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*J+l/2 


2 ^j+1 " ^ j+l/2^^j+l/2 


(3.8b) 


where 


-M 

'^ 1 + 1/2 


'^i+1/2 ^1+1/2 


^1+1/2 “ ®l^^^i+l/2 


3.8c) 


We can rewrite (3.8b) as 


^J+l/2 - 2 * h ‘Sj-^ «j+l- '’<^3+1/2-' V 1 / 2 ’ ''j+1/2 ''’• 


L«»TTmia 3.2 Suppose Q(x) is Lipschitz continuous and satisfies 
* ®j+l " “ ^''j+1/2^^^ •'^j+1/2 ^ ^ 


(3.9a) 


Tj+l/2*^j+l/2 ^ =- gjrl - = 0(A ) 


(3.9b) 


then the nunierical flux (3.8) satisfies (3.6) 

Pr oof: The modified numerical flux ^1|.3_/2 (-3*8d) differs from the original 

flux fj+Q^/2 ^ folla-ring way : 

^n/2 = * X '®j " Sj+i * 

(3.10a) 

“ ^^''j+1/2 ^ ^j+1/2^^ ^j+1/2 • 

Substituting (3.7) for fj+^/2 ^ (3.10a) we get that (3.6) holds if the 

relation 

^ * ^+1^ * ^^^'’j+1/2- ” ^^'’j+1/2 ■*■ ^j+1/2^^ ^j+1/2 ^ “ 

(3.10b) 

[Q(5.,1/2) - (v.,^/2>b V + O(i’) 


is satisfied. 
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Since Q(x) is Lipschitz continuous 




.< const. IY.+I/ 2 I, 


C3.10c) 


tiierefore if follows frcan (3.9b) that the second teim on the LHS of (3.10b) 

2 

is itself 0(A ); consequently (3.9a) iinplies (3.10b). This conpletes the proof 
of Lemm 3.2. 


We construct gj. = g(-V£_ 2 > tfiat satisfies (3.9) in the following 

way; 




Si-1/2 -Si+i/2)] 


(3.11a) 


0 


when » 0 

when gi.i/2‘Si-l/2 ° 


where 


Si+1/2 ■ 7 *-^^''i+l/2^ " ^'' 1 + 1 / 2 ^ ^ \+l/2 ^ 


(3.11b) 


®i+l/2 ” Cg,’4.T/9 5* 


i+1/2- 


(3.11c) 


Lemma 3.3 ; Let be defined by (3. 11) , then relations (3.9a) and (3.9b) are 
satisfied, and 


^^j+1/2^ ■ ^®i+l ~ ^ 7I ^^^^3+1/2^ " ^'’3+1/2^ 


(.3.12) 



- 17 - 


Proof ; First let us assume tha*- Pj+i/2*^_l/2 ^ using the definition 

<3. 11a) and the relation min (a,b) = j [(a+b) - |a-b|] 


■ ? ^®j-l/2 * Sj+i/2 “ “j+i/2^®j+l/2 “ ^;)-l/2l^ 

” ^±1/2 * ^’^®j+l/2 “ ®j-l/2^ " ^3+l/2'^j+l/2 " ^j-l/2l 


(3.13a) 


From (3.11b) we conclude that if v is smooth and Q(x) at least Lipschits 
continuous then 


^ ^ nfA^\ 

^+ 1/2 ’ ^- 1/2 - 


Thus (3.13a) and (3.13b) inply that 


gj = Sj±l/2 * 0^'' >• 


(3.13c) 


It is easy to see that C3.13c) holds even if Sj»i/2 ’^j+1/2 ^ then 

j = 0 but gj+i /2 = O(^i^) itself (since ‘^j+i /2 = 0(A^)). 

Relations (3.9a) - (3.9b) follow ijmediately by rewriting (3.13c) as 


gj = gj+i/2 O(A^) , = Sj+1/2 


We turn nw to Dix)ve (3.12). We observ'e from the definition (3.11a) that 


gj and gj+^ cannot be of different sign, hence 
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|gj+l - 8jl S n>»<<|gj|.|gj+i|) 

‘ iVl/2'- 

Thus it follows from (3.11b) that 

Itj.l/jl = Igj.i - gjl-'Ujn/j ''I ‘jn/2 ''I 

= IQ<*j+l/2l ■ ^'’j4l/2’^l ’ 

this ccfipletes the proof of C3.12). 

We shew now that the 5 -point second order accurate scheme (3.9) with (3.11) 
is TVNI under the same CFL restriction of the original 3-point first order 
accurate scheme (3.1). 

Lemma 3.t^; Suppose Q(x) satisfies (3.'?) and is defined by (3.11), then 
the finite-difference scheme (1.4) with the numerical flux (3.8) is TVNI under 
the CFL restriction (3.3). 

Proof : Since (3.8) is (3.1) applied to a modified flux (3.8a), it can 

— -M 

be rewritten as (3.5) replaced by '^j+1/2 ^3. 8b). Wa conclude 

from Lemma 3.1 that the scheme (3.8) is TVNI under the modified CFL 
restriction 
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, (3.14) 

nBXj|vV^l/2l J y 

To corplete the proof of Leinna 3.4 we shew that (3.14) is implied by the original 
CTL condition (3.3). Using (3.12) and (3.2) we get 


I ''j +1/2 1 ■ l'^j+1/2 * ^3+1/2! ^ Nj+i/2l l^j+1/2^ ^ ^''j+1/2* + 

+ I l^^''j+l/2^ ‘ ^'’j+1/2^^1 ^ l''j+l/2l ■*■7 ^ ■ 

= 1 ’ 7 ^l"j+l/2l “ 1 

vdienever l^j+i/2l ^ ^^^j+1/2^ ^ this shows that (3.. 14) holds. 


Reirarks ; (1) If Aj+ 1 / 2 ''' " ^ follows ijiinediately from the definition (3.11) 

that E. = g.j, = 0. This shews triat the modified nuineric-1 flux (3.Sd) is consistent 
*=3 “ 3+1 

with the physical flux f(u) in the sense of (1,4c). The schere (3.8) + (3.11) is TVin 
and • therefore it has convergent subsequences for all initial data of bounded 
total variation: the limits of these subsequences are week solutions of the scalar 
conservation law (2.1). To complete the convergence preof one has to show that 
all thiese limits are the same. In the constant coefficient case the solution to 
the IVP (2.1) is xnique and therefore the scheme is convergent (note tJat 
the scheme is nonlinear even in the constant coefficient case!) In the nonlinear 
case) com^ergence will fellow if one shows that the scher ? is consistent with 
Oleinik's entropy condition in the sense of (1.5). We shall discuss consistency 
vrith the entivpy condition in Section 5. 
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(2) Condition (3.6) is only a necessary condition for second order accuracy. 

It becomes a sufficient condition if the coefficient in the G(A^) term in (3.6) 
is differentiable, except possibly at a firi.te r.uirber of points N(t, At), such 
that NAt -^0 as At 0 for all t. It is clear f’xjm (3.13) and (3.9) that 
the troublesome points where the scheme (3.8) + (3.11) nay degenerate locally to 
OCa ) truncation error are those where is discontinuous, i.e. 

vhere Q(v) - v = 0 orUj^=0. The fact that the scheme is TVNI controls 
the possible increase of the number of local extremum points in the conputed 
solution. The schemes that we consider in section 5 all have the monotonicity 
property (see section 2); i.e. the number of local extremum points in the conputed 
solution is nonincr asing in time, and thus bounded by that of the initial data. 

(3) The modified equation of the scheme (3.1), i.e., the equation which it approxi- 
nates to second order accuracy, is (2.7) with 


A(u,X) = QCv) 


1 


V = Aa(u) 


(3.15a) 


We rewrite the modified ecuation as 


Ut + (f-^[Q(v) - v^] u^}j^ = 0, 
2 X 


(3.15b) 


and observe that g£+]^/2 ^3. 11b), and consequently (3.11a), is an 

approximation to the term 


g w ^X ” 7 tQ(v) - v^3 (Ax u^). 


(3.15c) 
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Our rrethod to convert a first order accurate IVNI scheme into a second order 
accurate TVIII scheme is based on the follcwing heuristic argument: The first order 

scheme C3.1) approjcimates 

+ Cf-^]„ = 0 (3.16a) 

X A .X 

to second order accuracy. Therefore, applying the same scheme to 

+ [f + ^3^ = 0 (3.16b) 

results in a second order accurate approximation to u^ + f^ = 0. To bs able 
to apply the scheme to (3.16b) we have to define g(u) with a bounded derivative 
dg/du, therefore g^+3^/2 ^^.llb) is replaced by (3.11a) (see C 5 3). 
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4. Systems of Conservation Laws . 

In this secticn ve describe how to extend our new scalar scheme of s^ion 3 
to systems of conservation laws. Our extension technique is a somewhat generalized 
version of the procedure suggested by P. Roe 5n [14]. The basic idea is to extend 
the scalar scheme to the system case by applyiny it”scalarly” to each of the (approp- 
irLately linearized) characteristic variables. 


Let S(u) = (R^(u), r2(u),...,R^(u)) 


(4.1a) 


be a BBtrix, the colcmns of which are the ri^t eigenvectors of the Jacobian 
matrix ACu) (1.2a). Then 


S"^AS = A 


A. . a^(u)6 . . 
11 11 


(4.1b) 


1 O' TTI 1 

The rows L‘‘‘Cu), L^(u) , . . . ,L Cu) of S”‘^(u) constitute a conplete system of 

left eigenvectors of A(u) which is bi-orthonornal to the system of right eigen- 
vectors, i.e. 


L^R^ = 


(4.1c) 


In the constant coefficient case A(u) s A = const. 


u. + Au = 0, u(x,0) = <j>(x), 


-00 < X < 


(4.2) 


one defines characteristic variables w = (w ) by 


k ,k. 
w = L u. 


w = S 


(4.3) 
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It follows (4.1) that (4.2) decouples into m scalar characteristic 

equations, 1 j k ^ m 

wj + = 0, w^(x,0) = L%(x), -• < x < •. (4.4) 

This offers a natural way of extending a scalar scheme to a constant coefficient 
system of equations (4.2) by applying it "scalarly" to each of the m scalar 
characteristic equations (4.4). 

k 

The characteristic variables w in (4.3) can also be viewed as the con>* 
ponents of u in the coordinate system {R }, i.e. 

u = ^ (4.5) 

k=l 

We use this interpi'etaticn of characteristic variables to extend the scalar 
scheme to general nonlinear systems of conservation laws. 

Let v^+i /2 = ^'j+ 1 ^ be an average of v^ and v^^^, i.e. a sirooth 

function V(u,v) such that 

V(u,v) = V(v,u) (4.6a) 


V(u,u) = u 


(4.6b) 
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V 

and let “j +^/2 denote t!.e c»nponent of in the coordinate 


system ^*^^^^+ 1 / 2 ^^ 


^j+1/2^ " ^+1/2 

(4.7a) 

k k 

® j+1/2 " 4+1/2 ^j+ 1 / 2 '' ’ 

(4.7b) 


here we use the notation convention " ^^''j+l/2^ “ bCVCv^, v.,i)). 


Vfe new extend the sc^ar schsne (3.8) + (3.11) to general systems c' 
conservation laws as follows; 


Vj ' ''j "^^^+1/2 " ^-1/2^’ 


(4.8a) 


^j+1/2 " ^ 


2T jJi -j+1/2^^ ® j+1 " ^ ^ j+1/2 ^j+1/2^ j+1/2^ 


(4.8b) 


Ic Ic 

where '^j+ 1/2 “ ^''j+1/2^ 


H ■ s^^i^/jnaxCO, inin(|g^^.j_/2i »gi_i/2*Si^l/2^^> ®i+l/2 ‘ (4.8c) 


. ,iVk 1 '>'k 


♦ /~k 


®i+l/2 "7 ^ i+1/2^ “ ^ i+1/2^ ^ i+1/2 


C4.8d) 


k 

^in/2= 


fjc kv . ?< 

%+l “ ^i^'^“i+l/2 


when “i +]^/2 ^ 0 


k 

‘i+1/2 


(4.8e) 


0 


vdien 
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The second order accurate one-step Lax-Wendroff sclieme can be represented 


as 


rt" 


j+l/2 ■ 7 


Cf(v.) + f(v.., )] - iij I “-i+i /9 


'i+i'-* T\ 


k=l 


j+l/2' "j+l/2 1+1/2 


(4.9) 


5 7 CfCVj) + - 7 l^A^''j+i/2^ ^^^j+ 1 / 2 '' 


Lenina 4.1. Suppose {Q^(x)} are Lipschitz continuous, then (4.8) satisfies 


^jU/2 " ^+1/2 


(4.10) 


Proof: Rawrite C4.8b) as 


^j+1/2 


vliere 


k 

■^j+1/2 


_ --U-/ . V pk k 

" ^ j+l/2 3 +1/2^ j+l/2 

" Yx ^j+l ■ ^^^j+l/2^ " 


/ Ic V 2 •! jc \ 

^ j+l/2 ' ^ “j+1/2' 


- Ij- 


k X ^k^ k V 1 k 

j+l/2 ' ' j+l/2 ^ ^ ^ j+l/2 j+l/2 ’ 


(4.11a) 


(4.11b) 


and then use (3.10c) and conclude (3.9) from Leima 3.3. 


Vfe define the total variation TV(v) of the vector mesh function v to 
be 

• m . 

TV(v) = I I ^‘*•12) 

j=-a> }c=l 

vtiere “^+2/2 defined by (4.7), and shov;: 

Leirirra 4.2 Suppose Q (x) satisfies (3.2) for all k, and that A(u) = A = constant, 
then the scheme (4.3) is TVTNT under the C?L restriction 
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XtnaxjA^I ^ y = min ^ 1 (4.13) 

V 

Inhere u are the restrictions in (3.2). 


Rpoof; Because of the assunption A(u) = const) 0^^} , (L^) and {a^} 

V 

are all constant. Multiplying (4.8b) from the left by L , we obtain (3.8d) 
for the characteristic variable w^ (4.3); and '^^ 1/2 ^ becomes 

identical with (3.11). Thus by Leima 3.4 we conclude that under condition (4.3) 
the total variation of each of the characteristic variables is nonincreasing, 
and ■therefore the total variation (4.12) is nonincreasing as well. 


Corollary 4.3. The scheme (4.8) in "the constant coefficient case is 
convergent under "the res-triction (4.13) for all initial data of bounded total 
variation, and is second under accurate. 

We renerk ■that this corollary is not -trivial since the scheme is highly 
nonlinear even in the constant coefficient case. 

Our technique to extend scalar schemes to the system case does not require 
any particualr form of averaging VCu,v) (4.6)* Roe in P-5 3 is using a specific 
foim of averaging -that on top of being mathemtically pleasing, also enables the 
CCTiputational advan-tage of perfectly resolving stationary discontinuities. 

In [73 we shCT7 -that if the system of conservation laws (1.1) x»ssesses 
an en-tropy function (.1.3), -then it is symme-trizable , and there exists a mean value 
Jacobian A(u,v) such ■that 
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(i) f(v) - f(u) = A(u,v) (v-u) 


(4.14a) 


(ii) A(u,u) = A(u) (4.14b) 

(iii) A(u,v) has real eigenvalues {a^(u,v)}|^_^ and a conplete set of right 
eignevectors {I^(u,v)}JJ.^ . 


In the context of the scheire (4.8) Roe’s extension technique is expressed by 

taking ^j+i/2 ^ ^^*7) - (4.8) to be the eignevalues and the right 

eigenvectors of the mean value Jacobian A(v^, Vj^^) (4.14aj, respectively. 

• Ic 

Thus if Oj+3^/2 defined by (4.7a) 



I 


k=l 


a 


j+1/2 -^j+1/2 


(4.15a) 


then it follows from (4.14a) t'nat 


m 



k k 
“j+1/2 ^ j+1/2 



(4.15b) 


The lelation (4.15) makes the scheme (4.8) a more faithful extension of (3.8) 
in the sense -that (4.8) for m=l is identical with the scalar scheme (3.8). 


Vfe observe that if f(u) = f(v) in (4.14a) then v-u is a ri^t eigenvector 

or A(u,v) , ccrresocnding to a zero eignvalue. Hence in (4.15) a (u,v) = 0 for 
k ‘ k 

k ^ k^ and a (u,v) = 0, for some k^. It is easy to see that if Q °(0) = 0 

in (4.8) then the stationary discontinuity 


u(x,t) = (x) 



X < 0 
X > 0 


, f(u) = f(v) 
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is £l1so a stationary solution of (4.8). (see CIS]). 

In the case of the Euler equations of gasdynamics , where the flux f(u) 
is an hcxnogeneous function of u of degree 1, it is possible to express A(u,v) 
in C4.14) as 

ACu,v) = A(VCu,v)). (4.17) 

This relatively siirple function V(u,v) (see [15]) will be described in 
Section 7. 

Reniarks; (1) Note that we use Q^(x), thus allc^/inp; different functions (3.2) 
for different characteristic fields. As observed by P. Roe [14] the extension 
technique of this section Permits even the use of conp].etely different scalar 
schemes for different characteristic fields. 

(2) In most applications it my be advantagous to replace the term 

^ Cf(Vj) + f(Vj^^)] in C4„8b) by = f(V(Vj, Vj^^)). For sure this 

sinplifies the prxDgramming and reduces the CPU time, without altering the min 
properties of the scheme. A possible disadvantage of such a change is that the 
scheme my lose the property of perfect resolution v^f stationary discontinuities. 

(3) The particular definition (4.12) of total variation is motivated by the 
definition of Glim’' s functional in [ 3 ] . When applied to a piecewise-smooth 
solution u(x,t) of (1.1) 

lim TV(u) = / I |L*^(u)u |dx + 1 ? ia^x*)| 

1<?1 J ]^1 ^ 


(4.18) 
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k k 

where Xj are points of discontinuity, and a (x^) denotes the value of 
in C4.7) evaluated v;ith respect to Vj^ = u ((x. )_,t), v^^^ = u((x.)^,t). 


There is no reason to expect that the functional (4.18), and consequently/ 
C4.12), is generally nonincreasing with t. Based on Glijmi’s results T3 ] we 
do however believe that this functional (under certain conditions) is bounded in 
t. At this time we do not have estijiates of the possible increase in total 
variation in solutions of the scheme (4.8), and therefore cannot prove convergence 
in the nonlinear sS”stein case. 


- 30 - 


5. On the choice of 0(x). 

In section 3 we have presented the basic idea of our new scheme in the 
scalar case, which can be described algorithmically by: "Take a 3-point TVNI 
scheme (3.1), vdiich is first order accurate, and apply it to the modified flux 
viLue (3.8) to obtain a 5-point second order accurate TVNI scheme (3.11) + (3.8) 

We begin section 5 by considering different choices of Q(x) in the 
3-point scheme (3.1). In the end of this section we argue that most of the pro- 
perties of the 3-point point scheme (3.1) go over to the irodified 5-poirit scheme 
(3,8) +(3.11). 


In section 6 we discuss the system case. 


A natural cohice of Q(x) in (3.2) is Q(x) = [x[ , which corresponds to 
the least dissipative TVNI scheme of the fom (3.1). Ihe scheme (3.5) with 
Q(x) = |x| can oe re^'?ritten as 


v?^^ = V? - (v.^, ,o) A. /«v^ - (v. , /«)^ A 


n 


3+1/2' “j+1/2 


j-1/2' ‘"j-1/2' 


(5.1a) 


v^ere 


\T = min(v,0> = ~ (v - |v|), v‘*’= max(v,0) = y (v + jv|). 


(5.1b) 


This scheme (5.1) is a generalization of the well-known i^stream 
differencing sdieme of Co\jrant, Isaacson and Rees, and it is well investigated 
in the literature (see [5], [10] and the references cited there). We new give 

a brief review of its relevant properties. 
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Leiura S.l» Solutions of the upstream differencing schere (5.1), under the CFL 
r e s triction 

Xuax |a(u)| ^ 1, ^ = range of initial data, : (5.2) 

uC<> 

have the monotonicity property stated at the beginning of Section 2. 

Sketch of Proof; Let V7(x) be the piecewise-linear interpolant of v” , i.e. 
V7(x.) = , X. = jAx. 

Then (5.1), under the CFL restriction (5.2), iirplies that there exists 
a monotone ncndecreasing sequence { Xj) such that 

= V/(x.), I S. - X. I ^ 

Ihis shCTvs that has no more local exTremum points then v^, and that 

the value of a local ndniirum (raxirnurn) is nondecreasing (nonincreasing) , which 
is the assertion of Ijerma 5.1. 

V7e observe that the stationery 3 urr!p discontinuity (4.16), admissible or 
inadmissible, is also a satationary solution of the upstream differecing scheme 
(5.1). On one hand this is a desirable propei-^ as it djrplies good resolution of 
stationar:.’ shodoion the other h^x^d it indicates that the scheme many select 
nonphysical weak solutions that do not satisfy th * entropy condition (1.3b). This 
property.' is related to the fact that the viscosity term 9(u,X) (3.15a) vanishes 

for V = 0, and it is conmon to all schemes (3.1) with Q(C) = 0. ’ 
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To prevent entropy condition violation of this sort and make Q(x) 
smoother at the same time, we modify QCx) = |x| near x = 0 to be 


QCx) = 



for |x| $ 2e 

for |x| i 2c 


(5.3) 


with say e = 0.05. Csee Fig. 1) 


This change has increased theamour.t or numerical viscosity*' for |x| < 2e 
so that now B(u,A) > 0 for lv| < 1. p(u,X) finishes only for v = 1; 

which can be handled by taking v < 1, say y = 0.95, in the CFL resti’iction 
(3.3). Consequently the scheme (3.1) with Q(x) C5.3), cannot have perfect 
resolution of “die stationary shock (4.16). 


Another possible choice of Q(x) which corresponc? to a scherie that has 

already been investigated is Q(x) = x^ + ^ (see C ^ 3 r.nd [ s3). This scl'iema 

is the Lax-Wendroff scheme modified by the addition of the linear viscosity term 

4<v.^, - 2v. + V. ,). Since now Q(x) - 1x1 = (Ixl % 0, it is the require- 

8 3+1 3 3-1 ' ' 

ment QOOj .l in (3.2) which restricts the CFL condition (3.3) to 
y = /? /2« 0.866 < 1 Csee Fig. 1). 

Vte turn now to discuss the second order accurate TVNI scheme (3.8) + (3.11). 
Vfe note that the truncation error of this scheme is O(i^). except possibley where 
= o, indeperxient of the particular form of Q(x). Thu: the modified scheme 
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2 

renoves all 0(A ) errors due to nurnerical viscosity, except at the points where 
goes through zero; there some viscosity is needed to ensure monotonicity. 
Consequently we expect the new scheme to have high resolution of shocks, stationary 
or moving, almost independently of the particular choice of Q(x). 

Based on these considerations it se< 2 ms to us that (5.3) is the best choice 
of QCx) for (3.11) 

Ifext we deal with the question v.’hether properties of the scheme (3.1) 
go over to the modified scheme (3.8) + (3.11). We expect properties of the original 
scheme that can. be expressed as a relation between F(v^) and F(v^‘^^), where 
F is some functional that depends on mesh values only, to go over to the modified 
second order accurate scheme. We have seen that this is true for the TVNI pro- 
perty; v?e can also prove that for ()(x) = |x| the modified scheme has the mono- 
tonicity property, and that in thie stationary discontinuit>' case (4.16) it behaves 
the sane way as the original schene. 


To support our conjecture in general, we note thax the modified scheme is the 
original schene applied to a mc-dified flux f + ^ g (3.1Cb). Thus if g were 
strictly a function of u, our statemerit becomes trivially true. Hovjever g in 
(3.11), which is consistent with zero flax, iray be cosidered g(v^) only where 
is monotone. The particular form used in (3.11) is designed to match the 
different definitions of g(’/‘) in neighbour'ing monotone sections by setting 
g = 0 at local extremum points (see [ 5]). That this netching is smooth, is 


evident by the fact that y = Ag/du 


is w’ell-de fired eveinwhare, ana is bovffided in 


absolute value by (3.12). 
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Hopefully the piecewise definition of g(v^) does not spoil validity of 
such properties for the modified scheme. 

Based on these arguments we form the conjecture that if the original 
scheme C3.1) is consistent with the entropy inequality (1.3b), so is the modified 
second order accurate schane C3.8) + C3.ll). 
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6. On resoltJtion and entropy . 

The extension technique of section 4 consists of applying scalar schemes, 
not necessarily the same, to characteristic fields. Therefore we have the possibi' 
lity to custon fit the scheme to the computational needs of each characteristic 
field. 


Vfe consider here systems of conservation laws where the characteristic fields 
are either genuinely nonlinear (a|^R^ i 0) or linearly degenerate = 0» 

see [11]). The waves of a genuinely nonlinear field are either shocks or rarefac- 
tion waves, depending whether the characteristics are convergent or divergent . 

The waves of a linearly degenerate field are exclusively contact discontinuities. 

First let: us consider the latter' case of a linear'ly degenerate characteristic 
field. It is well known that the Rieirann invariants of this field, one of wliich 
is the characteristic speed, are continuous across a contact discontinuity. There- 
fore the propagation of a contact discontinioity is ccmputa-ionally equivalent 
to that of a scalar discont,inuity moving with a constant characteristic speed. 

As renarked in section 3, the solution of the first order accurate scheme (3.1) 
is the same, uo to second order terms, as the solution of the modified differen- 
tial equation <3. 15b) Cor (3. 16a)). Tlie modification of the flux f(u) = constant, u 
by the addition of the term C- ^ g) has the effect of making the characteristic 
field slightly divergent in a region of a discontinuity. Ct«nsequently the coirputed 
discontinuity is being spread at the rate of where n is the number of 

time-steps taken j similarly, a p-th order accurate standard scheme spreads a contact 
discontinuity at the rate of (see [5 J). 
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Our new second order accurate scheme is obtained by applying the original one 
to the flux f + Y g (3.16b). This has the effect of reducing the divergence of 
the characteristics induced by the original scheme in a region of discontinuity, 
and we e>!pect the corputed discount inuity to spread like n^^^, rather than n^^^. 


To ensure that the computed contact discontinuity does not. spread in time at 
all, we want to make the effective characteristic field sli^tly convergent. To 
accorplish that we increase the size of g in (3.16b) e.g. by multiplying the RHS 
of C4,8c) by a + 6^), i.e. 


= C 1 + ^ ^l®l+l/2l ’^i-l/2*®i+l/2^^ 


where 




(6.1a) 


(6.1b) 


k 

“i+1/2 


is defined by (4.7). 


In [6 J we show that 


0 s ej ^ 1 


(6.2a) 


= OCax) in regions of smoothness. 


(6.2b) 


and that 8^ = 0(1) in regions of discontinuities. 


It is easy to see that this change in (4.8c) modifies the numerical fliox 
2 

(4.8b) by 0(A ) cind thus does not spoil the second order accuracy of the scheme. 


From (6.2a) and(3.12) it follo;^^ that now 
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(6.3a) 


Hence for Q(x) (3.2) we get in Lemna 3.4 


I '’j +1/2 1 ■^lYj+ 1 / 2 ! ^ I +1/2 1 ■*■ Q^''j+l/2^ " ^''j+1/2^' 


(6.3b) 


For Q(x) = |x| we have that the RHS of (6.3b) is less or equal 1 for 
jv^+ 1 ^ 2 l ^ consequently the modification (6.1) does not alter the CFL 

condition for this field. The same is true for Q(x) in (5.3) with e< 1/3. 
(Kos-jever for Q(x) = x + 1/4 (6.3b) inplies a more stringent restriction, 

r'amely hhxI ^ 3/4). 


We rerark that there are no entropy considerations associated with a 
linearly degenerate characteristic field. 


We turn now to the case of a gen'iinely nonlinear characteristic field, 
Where the corrputational aspects to be considered are resolution of shocks and 
enforcement of the entropy condition. 

Unlike contact discontinuitieL . shocks are formed and main- 
tained by local convergence of the characteristic curves. The 
reduction in numerical viscosity due to the addition of (4.8c) 
is usually sufficient to ensure good resolution of shocks. When 
higher resolution is desired/ one may employ a mechanism of the 
form (6.1) to enhance the local convergence of the characteristics 
in the shock region and thus improve its resolution (see [5]). 
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The conputational aspect of the entropy condition may be demonstrated by 
considering a Rieirann problem for (1.1) involving two nearby states and Uj^. 
lax in [113 proves (assuming all fields to be genuinely nonlinear) that there 
exists a unique "physical" solution to this problem, vMch consists in general ■ 
of both shocks and rarefaction waves. One can show, using exactly the same method 
of proof, that there also exists a weak solution consisting only of shocks , which 
may differ from the physical one. Whenever the first cxjntains rarefaction waves, 
the has nonphysical "expansion shocks". Therefore to ensure that the scheme 

selects a physically relevant weak solution, consistency »ith the entropy condition, 
on top of consistency with (1.1), is required. 

In section 5 we have argued that the scalar scheme (3.8) + (3.11) is consistent 
vdth Oleinik's entropy condition. The question arises now whether the system 
version (4.8) of this scheme will always select the physical weak solution. 

T,P,Liu 1 123 shews that if a weak solution of (1.1) contains only "admissible" 
discontinuities that satisfy a "scalar" Oleinik conditiu: with respect to the 
Rankine-Kungeniot curve, tiien this solution also satisfi:-5 the entropy inequality 
(1.3b). Using "diis theory the above question can be refcrraulated as: Will the 

scheiae C4.8) reject discontinuities llTat do not satisfy Aleinik's condition along 
the Fankine-Hugoriot curve? It seems possible to answe? this question affirmatively 
far sufficiently weak shocks, and feasibly for a larger class of problems. 

1 We remark that there is also a one-parameter family rf states connecting u^^ 
end through rarefaction curves only (see H3J). ferever a "negative rare- 

faction wave" is not corputationally realizable as it ir a multi-valued function 


of X 
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7. Application to Euler equations of Fasdynamics 

In this sectioii we describe the application of our new scheme (4.8) to 
the Euler equation of gasdynamics: 


+ f(w)j^ = 0 


(7.1a) 



■p 


'0 ’ 

w = 

m 

IE 

, f (w) = uw + 

P 

Ipw 


(7.1b) 


p = (y - 1) (E - 1/2PU ). 


(7.1c) 


Here p, u, p and E are the density, velocity, pressure and total energy, 
respectively; m = pu is the momentum and we take y = 1.4. 


The eigenvalues of the Jacobian matrix A(w) = f^ are 


a^Cw) = u - c , 32 (w) = u , 33(w) = u + c 


(7.2a) 


1/p 

where c is the sound speed, c = C Yp/p) 


The corresponding right eigenvectors are 



■ 1 


' 1 ■ 


■ 1 

R.Cw) = 
1 

u - c 
H - Vic 

4 4 

, R2(w) = 

u 

L 2 J 

, = 

u + c 
H + uc 


(7.2b) 


MO - 


i^ere H = (E + p)/p = c^/(y-1) + 

Let a^Cvj^ ^ ” ^»2»3> 

of linear equations CM. 7) 

Wj^ - Wl = o^^(V(Wj^,Wj^)) 


1 . 2 

^ u is the enthalpy. 

^ • 

be tJie solution of the following system 


(7.3a) 


vdiere V(wj^, Wj^) (M.6) is some average state; denote its velocity and sound 

V 

speed by u and c, respectively. To calculate a in (7.3a) we first 
evaluate 


= (y-1) {[E] + ^ u^Cp] - u [n0}/c2 


(7.3b) 


= {[m] - u[p]>/c 


(7.3c) 


vSiere [bj denotes [b] = - bj^; then a ' in (7.3a) are obtained by 


1 - I 


= 7 CC^ - C2), = [p] - C^, a = 7 (C^ + C2). 


(7.3d) 


The second characteristic field corresponding to the eigenvalue u is 

2 2 

linearly degenerate, i.e. , a^ R = 0; The other characteristic fields corres- 
ponding to the eigenvalues u ± c are genuinely nonlinear. Therefore in 

2 13 

(6.1a) we define ®j+3^/2 (6.1b) for k = 2 and set ®j+i/2 ~ ®j+l/2 “ 


Given w? we now list the operations needed to calculate 


0 . 
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(i) Cbnpute: (Uj = rij/Pj) ^lnd (Pj) by (7.1c) 

(ii) Coffipute ^j+1/2 calculate 

■*■ °j+l/2^’ evaluate k = 1>2,3 by (7.3b) - 7.3d). 

Define X = At/Ax = v/maXj (| ■•■ ^j+i/2^ where \i is the prescribed 

CFL restriction in (4.13) 


(iii) Conroute: ~ ^^j+1/2 (7.2a); Sj+2/2 (^.8c) 


®j+l/2 (B-lt>) for k = 2; Set 6j+3^/2 = ^jn/2 


= 0 . 


(iv) Calculate by (6.1a). 


(v) Corrpute Yj+1/2 .8e), and by (4.Sb) and relations (7.2b) 


(vi) Confute by (4.8a). 


We note that in corq^uting ^j+3_/9 ©‘"i- ccvild ta]-:e advantage of the 

V 

sinple form of the R in (7.2b). 


Ne>rt: how to inpleinent Roe’s linearization technique (4.14) - 

C4as) in the above algorithm. Roe presents a pai-ticular form of averaging 
V(»L ,v/^) such that for the Euler equations of gasdyiiamics , the nean value 
Jacobian ACw. ,w„) in (4.14a) cai be e:<pressed bv (4.17). Tliis averaging takes 

ti XV 

the followduig fom; 
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Vl/2 . 

(7.4a) 

®j+l/2 ■ ^^Y-1) (flj+i/2 - 
\ihere <b> denotes the arithemetic mean 

<'b> = 7 (bj + b^^3^) (7.4b) 

Therefore to use Roe’s linearization in our scheme all one has to do is to 
ccnpute and ^ (7.4). 

Vfe renark that the averaging in C7.4) is rather expensive. It seems to us 
that in nany applications the sinple arithmethic average (7.4b) vd.ll do just as 


vfell. 
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8. !?unierical experijnents . 

In this section we present sane numerical experiments that demonstrate 
the nerfomance of the proposed second order accurate scheme. Vte consider here 
the following two versions of it: 



Vj - 



) 


^j+1/2 = 2 + f(Vj+]^) - " gj+1/2 Rj+i/2^ 

where 

®j+l/2 ■ ^^'Vl/2 * +1/2 +1/2 “ ^ ^j+1^ 


^8. la) 
(8.1h) 
( 8 . 2 ) 


and ^ is defined by an appropriate variant of (6.1) and (4.8d); 
y^^2./2 defined by (4.8e), 

The first version is (8.1) - (8.2) with 


C^bc) = lx] C8.3a) 

and will be referred to as the scheme ULTl; The second version is (8.1) - 
C8,2) with 

Q^(x) = x^ + 1/4 (8.3b) 


and will be referred to as the schene ULT2. 

For conparison sake ve also present calculations with the folla'/ing two 
scheres : 
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(i) The second order accurate Lax-Wendroff-type scheme (8.1) with 
vAiich is referred to as the LW scheme. 

(ii) The first order accurate Godunov-type scheme of Roe (see [15] and [10]) 
%rfdch is defined by (8.1) with 

®j+l/2 = I ''j +1/2 1 

and is referred to as the ROE scheme. 


In all the schemes and experiments reported herein we use the Roe 
linearization (7.3) - (7.4). 


(I) The shock tube problem . 

Vfe cosider now a Riernann problem 


w(x,0)=J 


"l 

X 

o 

V 

"r 

X 

> 0 


(8.6a) 


for the Euler equations of a polytropic gas (7.1). Our first set of data is 



’ 0.445 


0.5 


0.3111 
.8.928 , 

II 

0 

1.4275 

4 / 


(8.6b) 
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Other nunerical experiments with this problem are reported in [6 ] and the 
references cited there. 

In figures 2a, 2b and 2c we show the results obtained by the ROE, LW and 
ULTl schemes, respectively. The numerical values are shown by circles; the 
e^ct solution is shown by the solid line. The calculations in figure 2 were 
performed with 100 time-steps under the CFL restriction p=0.95 in (4.13), 
and 140 cells. 

In figure 3 we repeat the calculations presented in figure 2 for a diffeiont 
set of data for the Rieirann Problem (8.6a) 



’ 1. ■ 



0.125' 

h- 

0. 

2.5 

k 4 

» 


0. 

0.25 

V # 


Other numerical experiments with thi.s problem are presented in [i®]. The 
calculations in figure 3 were performed with 50 time-st-=ps under the CFL restriction 
y = 0.95 in (4.13), with 100 cells. 

\4e renark tlat the solution of ULT2 for these probelms look almost identical 
to those of ULTl. In both schemes we find almost no dependence on the CFL 
number. 


(II) The Quasi 1-D nozzle oroblem . 

Vfe consider an axis\TTmetric nozzle with a cross-section area A(x). The 
cross-section average of the flow satisfy the following one-dimensional system 
of equations 
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= - s(V7,x) 


s(w,x) = 


0 

dA 
P 3x 


1 


where w, f(w) and p are given in (7.1). 

2 

In figures 4,5 cuid 6 we present numerical approxinations to steady state 
solutions of (8.7a). 


In figures 4 and 5 we show solutions for a diver'gent nozzle vjith the cross- 

s^tion area 


A(x) = 1.398 + 0.347 tanh (0.8 m - 4); 


(8.7b) 


the flow condition is supersonic at the entrance and subscni at rhe exit. 

Figures 4a and 4b sho^‘7 steady state solutions on a cri.:de ?nesh of the ROE and 
ULTl schemes, respectively. Figure 5 sliows the ULTl results for the same problem 
on a finer mesh. 


In figure 6 we shew a steady state solution of *he ULTl scheme for a 
oonvergent-<iivergent nozzle with the cross-section area 


ACx) 


1 + (Aq - 1 ) (1 - x/5)^ 

1 + (Aj, - 1) [Cx - 5)/(x^ - 5)]^ 


X 5 
X > 5 


(8.7c) 


where A^ = entrance area, A^ = 
entrance as well as at the exit. 


exit area; 


Here the flew is subsonic at the 


2 These figtires are by courtesy of Helen C. Yee of the KASA-Ames Research Center. 
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The exact solutiorBin figures 4 to 6 are shown by the solid curve; the 
values of the numerical solutions are indicated by a rcrt)UP. 

Hopefully a detailed report of these calculations and the particular 
approach to steady state will be published elsev^re. 

(Ill) 2> D Flow throu^ a duct . 

In figures 7a and 7b we show solutions to the problem of the flow of air 
throu^ a duct contadning a step. Initially the flow is everywhere to the right 
at Mach 3, with p = 1.4, p = 1 and c = 1. The duct v;idth is 1, its length'is 3, 
and the step of hieght 0.2 is located a distance of 0.6 from the entrance. 

Figure 7 shows the results at t = 4 with a crude unifonn Cartesian grid with 
AX =Ay = 0.1. 

This problem was used by V/oodward and Colella to test the performance of various 
numerical schemes (see [IS] and the references cited htere). 


The solutions in figure 7 were obtained by a Strang -type dimensional 
splitting of the fom 





C8.8a) 


^ H- Hr He 


(8.Bb) 


where and 
approxuivating 



are one-dimenional finite difference operaters 


- 48 - 


♦ f(w)jj = 0 , Ly: + gCw)y = 0. (8.8c) 

If amd ly are stable and dissipative second order accurate approximations 
to the one-dijnensional equations in (8.8c), then the scheme (8.8a) - (8.8b) is 
a stable second order accurate approxijiation to the 2-diaensional problem. 


t f(w)jj + gCw)y = 0. 


(8.8d) 


In figure 7a we show for caiparison sake the results of the second order 
accurate hybrid scheme (see [6]) (8.1) with 


-k _ Ff k v2 , 1 ftk -* k 

®j+l/2 “ ^^'’j+1/2^ 4 ®j+l/2-‘ “j+1/2 


(8.9) 


vhere 9j+i/2 (6,1b) for all k. 


In figure 7b we show the results of ULT2. 


Both figures 7a and 7b shew 30 equally spaced density conto’jrs; Both 
calculations were perfonned with a CFL restriction of C.75 
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We rerark that the comer of the step in our calculations is treated as 
a shai"p comer without any rounding (or equivalent addition of numerical viscosity). 
The sonic line emerging from this comer is a curve on which an eigenvalue of (8.8d) 
vanishes. It is interesting to note that the results of the Godunov scheme in 
figure la of [18] indicate that part of the sonic line may turn into an expansion 
shock. We find a similar behaviour in the results of ULTl for this problem. 

Altogether v;e find the performance of the new second order accurate scheme to 
be quite pleasirig. Vte note that the scheme is siirple to pregram and requires only 
slightly more CPU time than a Lax-Wendroff scheme with some artificial viscosity. 
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Figure 2a . ROE scheme for (8.6b) 
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Figure 2b . LW scher.e for (8.6b) 
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Figure 2c. ULTl sdiene for (8.6b) 
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Gchene for (8.7b). 
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rinure 6. ULTl schcane for (8.7c). 
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Fig. 7b. ULT2 scheme - 10 x 30 grid. 
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